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Introduction 


Normalized Snow BRDF Models 


Comparison Summary 


Satellite remote sensing of clouds requires an accurate estimate of the clear-sky radiances for a given 
scene to detect clouds and aerosols and to retrieve their microphysical properties. Knowing the spatial 
and angular variability of clear-sky albedo is essential for predicting clear-sky radiance at solar 
wavelengths. The Clouds and the Earth's Radiant Energy System (CERES) Project uses the near- 
infrared (NIR; 1 .24, 1 .6 or 2.1 3 pm), visible (VIS; 0.63 pm) and vegetation (VEG; 0.86 pm) channels 
available on the Terra and Aqua Moderate Resolution Imaging Spectroradiometer (MODIS) to help 
identify clouds and retrieve their properties in both snow-free and snow-covered conditions. Thus, it is 
critical to have reliable distributions of clear-sky albedo for all of these channels. In CERES Edition 4 
(Ed4), the 1.24-pm channel is used to retrieve cloud optical depth over snow/ice-covered surfaces. 
Thus, it is especially critical to accurately predict the 1 ,24-pm clear-sky albedo a and reflectance p for a 
given location and time. Snow albedo and reflectance patterns are very complex due to surface texture, 
particle shapes and sizes, melt water, and vegetation protrusions from the snow surface. To minimize 
those effects, this study focuses on the permanent snow cover of Antarctica where vegetation is absent 
and melt water is minimal. Clear-sky albedos are determined as a function of solar zenith angle (SZA) 
from observations over all scenes determined to be cloud-free to produce a normalized directional 
albedo model (DRM). The DRM is used to develop a(SZA=0° ) on 10’ grid for each season. These 
values provide the basis for predicting p at any location and set of viewing & illumination conditions. 
This paper examines the accuracy of this approach for two theoretical snow surface reflectance 
models. _ , . 

Objectives 

• Determine errors in CERES Ed4 predicted clear-sky 1 .24 pm reflectances over Antarctica 

Data 

• CERES MODIS products derived using Ed4 algorithms: 

Aqua: 9-11 January 2008; Terra: 3, 15, and 20 January 2008 


A pproach 

Clear-sky reflectance 

P ~ o’ M-> — Xm(\^o’ 40 ^(Po)> 


( 1 ) 


p 0 = cos(SZA), p = cosine(viewing zenith angle VZA), <j> 0 = relative azimuth angle, x is a normalized 
bidirectional reflectance distribution function (BRDF), a = SZA-dependent clear-sky albedo. 

DRM 

S{[i 0 )= «(p 0 ) / «(p o =0° ) (2) 

Overhead-sun albedo is estimated for a given area by measuring the clear-sky reflectance, then 
computing the albedo using (1), and applying the DRM from (2) to solve for a(p o =0° ). Empirical 
(global averages) and theoretical (RTM) quite different (see below). 

BRDF Models: ^ m (p 0 , p, <|>) = p m (p 0 , p, <|>) / a m (p 0 ) 

p m { p 0 , p, <)>) computed using a radiative transfer model ( RTM, adding-doubling or DISORT) with the 
optical properties corresponding to snow-free observations. Snow surface is assumed to be a layer 
of randomly oriented ice crystals with optical depth of 1000. Model accounts for atmospheric and 
surface absorption to yield a TOA reflectance. 

Model A (used in Ed4)\ Single hexagonal column: length/width ratio = 750pm/160pm, effective 
particle diameter, D e ~ 310 pm ( Trepte etal. 2002). Regional (10’) values of a(A,,cj);p o =0 o ) updated 
every 2nd day to keep up with observed changes in the clear-sky reflectance fields. X = lat, (j> = Ion. 

Model B: Mixed habits of hexagonal column and plate, bullet rosette, and aggregate, which vary 
with particle maximum dimension D m . Surface roughness is 0.1. D e ~ 200 pm for the mixture. 

Estimating p m . / * Mo * l-M’o) = ,P 0 ) ^(l-O Xm( M-o’ 1-9 4*) 

Directional Models 
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Model A produces more forward and side scatter than Model B, which generates more backscatter. 

Overall, Model B has the more complex scattering patterns. 

Overhead-sun Clear-sky Albedos over Antarctica, 2008 


Model A Start Up 
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Model B Start Up 
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Model A updating tends increase the overhead albedo everywhere. Model A values are larger because the 
model-A DRM is nearly flat. The lower model-B values are much lower because the DRM increases with SZA. 


Results 

Instantaneous Comparison Example 


Terra 
1305 UTC 
9 January 
2008 
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Model A yields much lower reflectances in center of image and higher values near the ends of the 
granule. Model B looks more like the observations except in lower left corner. 


For model A, a large portion of 
the domain has differences 
between ±0.04, with more 
areas on the positive side. 
Extreme differences are seen in 
eastern Antarctica north of 
75° S. Overall positive bias. 

For model B, more areas having 
differences ±0.04 with a better 
balance than model A. Error 
pattern is substantially different 
from that of model A. Many of 
the largest errors occur in areas 
with significant altitude gradients 
or mountain ranges. 


Mean Reflectance Difference 
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Reflectance Difference (Model A - Observations) Reflectance Difference (Model B- Observations) 
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vza OO-K) 0.027 ± 0.045 
vza 10-20 0.023 ± 0.045 
v/a 20-30 0.017 ± 0.048 
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DRM accounts for SZA dependence; BRDF accounts for VZA and RAA 
variation. Errors in the predicted reflectance clearly vary with VZA. For 
model A, the bias increases monotonically with VZA from -10% at nadir to 
+1 9% at VZA = 55° . Standard deviation (SD) of the differences rises from 
1 1 to 1 9%. For VZA > 50° , a large portion of the differences 0.2 or 40%. 

Model B errors are less extreme with biases dropping from 6% at nadir to 
0% at 55° . The SD values are relatively steady with VZA at ~1 1 %. Few 
errors > 0.2 are seen for this model. 

Larger model A biases suggest DRM is not accurate. Error variations with 
RAA and SZA need further analysis. 


Ratio of Obs/Pred Difference of Obs - Pred 



MEAN 

STDEV 

MEAN 

STDEV 

Model A Predicted 

1.1099 

0.5639 

0.0442 

0.0988 


Model B Predicted 1.0541 1.1897 0.0067 0.0591 


Summary & Future Work 

• Model B clearly represents reflectance field better than model A. Model A optimizes mid VZA regions. 

• Future work will focus on adjusting model B based on Aqua and Terra observed clear snow 
reflectances to minimize SZA, VZA, and RAA-dependent errors 

- errors to be assessed for other areas: Arctic sea ice and land, IGBP types with seasonal snow 

- sastrugi and terrain roughness will also be considered in model adjustment 

• Revised model will be used in future CERES Editions 
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